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ISOTHERMIC TRIANGULATED SURFACES 


WAI YEUNG LAM AND ULRICH PINKALL 


Abstract. We found a class of triangulated surfaces in Euclidean space which 
have similar properties as isothermic surfaces in Differential Geometry. We 
call a surface isothermic if it admits an infinitesimal isometric deformation 
preserving the mean curvature integrand locally. We show that this class is 
Mobius invariant. Isothermic triangulated surfaces can be characterized either 
in terms of circle patterns or based on conformal equivalence of triangle meshes. 
This definition generalizes isothermic quadrilateral meshes. 

A consequence is a discrete analog of minimal surfaces. Here the Weier- 
strass data needed to construct a discrete minimal surface consist of a trian¬ 
gulated plane domain and a discrete harmonic function. 


1. Introduction 

Isothermic surfaces are central objects in classical Differential Geometry. They 
include all surfaces of revolution, quadrics, constant mean curvature surfaces and 
many other interesting surfaces m- In particular, all classes of surfaces that are 
describable in terms of integrable systems in some way or other seem to be related 
to isothermic surfaces [121 US] ■ 

A smooth surface in Euclidean space is called isothermic if it admits confor¬ 
mal curvature line parametrization around every point. Note however that there 
are various characterizations of isothermic surfaces that do not refer to special 
parametrizations. 

Discrete Differential Geometry lies between Discrete Geometry and Differential 
Geometry. The geometry of a discrete surface is determined by the positions of a 
finite number of vertices, such as those of a triangulated surface in Euclidean space. 
Smooth surfaces in Differential Geometry can be regarded as limits of discrete 
surfaces by refinement. The goal of Discrete Differential Geometry is to look for 
mathematical structures on discrete surfaces as rich as their smooth counterparts. 
It has many applications, for example in computer graphics and architectural design 

m- 

The same spirit applied to complex analysis has led to two different definitions 
of conformality for planar triangular meshes. One of these two is the theory of 
circle patterns [36] . where the conformal structure is defined by the intersection 
angles of neighboring circumcircles. It is motivated by Thurston’s circle packings 
as a discrete analog of holomorphic functions [34]. Another version of discrete 
conformality is based on conformal equivalence of triangle meshes |29[ 139] , where 
the conformal structure is defined by the length cross ratios of neighboring triangles. 
Luo introduced this notion when studying a discrete Yamabe flow. Its relation to 
ideal hyperbolic polyhedra was investigated in [5|. 

Previous definitions of discrete isothermic surfaces were all based on quadrilateral 
meshes that provide a discrete version of conformal curvature line parametrizations 
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of isothermic surface Emu]. Inspired by discrete integrable systems [5], Bobenko 
and Pinkall [2] considered quadrilateral meshes with factorized real cross ratios, 
which led to further investigation of discrete minimal surfaces and constant mean 
curvature surfaces [50] . Recently, the notion of curvature was introduced to discrete 
surfaces with vertex normals 13121]. 

Here we aim for a definition of isothermic triangulated surfaces which does not 
involve conformal curvature line parametrizations. It is motivated by a known 
(although not well-known) characterization that a smooth surface in Euclidean 
space is isothermic if and only if locally it admits a nontrivial infinitesimal isometric 
deformation preserving the mean curvature. The only reference that we could find 
is from Ciesliiiski et al. m, stating that this theorem was known in the 19th 
century. 

Infinitesimal isometric deformations of triangulated surfaces have been exten¬ 
sively studied since Cauchy’s rigidity theorem of convex polyhedral surfaces min]. 
An infinitesimal deformation of a triangulated surface in space is an assignment 
of velocity vectors to all the vertices. We can then calculate the change of edge 
lengths. An infinitesimal deformation is called isometric if the edge lengths are 
preserved. 

Suppose we have a realization / : H —>■ of a triangulated surface M = 

(V, E, F) such that each face of / spans an affine plane. Given an infinitesimal iso¬ 
metric deformation / : H —>■ R^, each triangular face {ijk} rotates with an angular 
velocity given by a certain vector G R^. These vectors satisfy a compatibility 
condition on every interior edge {ij}'. 

(1) df(^Cij') — df(^Cij^ X Zjjk — ^ ^jilj 

where {ijk} G E is the left face of eij and {jil} G F is the right face. 

On the other hand, it is well-known that the integral f F[ dA of the mean curva¬ 
ture has a very canonical discrete analogue "Yf Hy . Here we have defined the mean 
curvature associated to edge {ij} as 

Flij := (y.y\df{eij)\ 

where is the dihedral angle at the edge {ij} [40]. Under the infinitesimal iso¬ 
metric deformation given by Z on faces (Equation Q), we have 

Hij — d^ |c(/(e^)| — (^df (^eij{^ Zy]^ Zjiij. 

If we further demanded Fly = 0 on every edge {ij} then the infinitesimal isometric 
deformation would be trivial, i.e. an infinitesimal Euclidean deformation. Hence 
we consider instead the change of the integrated mean curvature around vertices 

Hi '■= ^ , ^ij\df {ey)\ = '^^{df (cy), Zyk — Zjii). 

3 3 

We are now ready to define isothermic triangulated surfaces. The smooth coun¬ 
terpart of the following formulation for isothermic surfaces is given by Smyth |38j . 

Definition 1.1. A non-degenerate realization / : U — )■ R^ of an oriented trian¬ 
gulated surface, with or without boundary, is called isothermic if there exists a 


R^-valued dual 1-form 

r : —)• R^, not identically zero, such that 

(2) 

^T(e;-) = 0 V^GU„t 

(3) 

J 

df{ey) X T(e* ) = 0 V{*j} G Ei^ 

(4) 

'^{df{ey),T{el)) = 0 Vf G Unt 
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Here and Vint denote the set of interior oriented dual edges and the set of 
interior vertices of M. 


The following is an immediate consequence of our definition. 

Corollary 1.2. A strongly non-degenerate realization of a simply connected trian¬ 
gulated surface is isothermic if and only if there exists an infinitesimal isometric 
deformation that preserves the integrated mean curvature around vertices but is not 
induced from Euclidean transformations. 

We will state several results about isothermic triangulated surfaces that closely 
reflect known theorems from the smooth theory. In Section and we prove 

Theorem 1.3. The class of isothermic triangulated surfaces is Mobius invariant. 

Theorem 1.4. For a non-degenerate realization f : V ^ of a closed genus-g 
triangulated surface the space of infinitesimal conformal deformations is of dimen¬ 
sion greater or equal to \V\ — 6g -\- 6. The inequality is strict if and only if f is 
isothermic. 


Theorem 1.5. Suppose f : V ^ is a non-degenerate realization of a simply 
connected triangulated surface. Then f is isothermic if and only if there exists 
an infinitesimal deformation that preserves the intersection angles of neighboring 
circumcircles and neighboring circumspheres but is not induced from Mobius trans¬ 
formations. 


Note that Theorem |1.4| concerns the theory of conformal equivalence of triangle 
meshes [5S1[33] while Theorem 1.5 deals with the notion of circle patterns [55] . 

In Section [^ we show that our definition generalizes isothermic quadrilateral 
surfaces |5]: Subdividing any isothermic quadrilateral surface in an arbitrary way 
we obtain an isothermic triangulated surface. 

In Sections mi and we provide examples of isothermic triangulated surfaces 
that are not obtained via quadrilateral isothermic surfaces. Triangulated cylin¬ 
ders generated by discrete groups as well as certain planar triangular meshes and 
triangulated surfaces inscribed in a sphere are isothermic. 

In Section[^we introduce discrete minimal surfaces via a discrete analogue of the 
Christoffel duality. Our discrete minimal surfaces are obtained as the reciprocal- 
parallel meshes for triangulated surfaces (with boundary) inscribed in the unit 
sphere. This approach mirrors the property that a smooth minimal surface is 
a Christoffel dual of its Gaufi map. The Weierstrass data needed to construct a 
discrete minimal surface consist of a planar triangular mesh and a discrete harmonic 
function. Such harmonic functions were first introduced by NcNeal [30] . They were 
used in linear discrete complex analysis since DufSn |16j and have applications in 
statistical mechanics (see Smirnov m)- 

In Section El we review the smooth theory and prove some new theorems that 
are similar to discrete results established in earlier sections. 

Throughout we use the language of discrete differential forms and quaternionic 
analysis as introduced by Desbrun et al. m and Pedit and Pinkall [31] . 


2. Notations 

Definition 2.1. A triangulated surface M = (V, E, F) is a finite simplicial complex 
whose underlying topological space is a connected 2-manifold with boundary. The 
set of vertices (0-cells), edges (1-cells) and triangles (2-cells) are denoted as V, E 
and E. 

Without further notice we assume that all triangulated surfaces under consider¬ 
ation are oriented. 
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Definition 2.2. A non-degenerate realization of a triangulated surface M in is 
a map f : V which is linear on each face and fi ^ fj for every edge {ij} G E. 

We say / is strongly non-degenerate if every face of / spans an affine 2-plane. 

We denote V^t and Ei^t the set of interior vertices and the set of interior edges 
respectively. We write Cy as the oriented edge from the vertex i to the vertex j. 
Note that Cy ^ Oji. The set of oriented edges is denoted by E. The set of interior 
oriented edges is indicated by Eint- 

We recall some notions about discrete differential forms [T5]. A (primal) 1-form 
uj : E is a. function defined on oriented edges of M such that 

A 1-form (jj is closed if for every face {ijk} € E 

uj{ey ) -\- uj{ejk) + i^{eki) = 0 . 
ft is exact if there exists f : V such that 

■“ fj fi ~ ^i^ij)- 

ft is easy to check that exactness implies closedness while the converse holds if the 
discrete surface is simply connected. 

Similarly we consider a 1-form r : E*^^ —> K on the dual cell decomposition 
M* = {V*, E* ^ F*) of M and call t a dual 1-form on M. Here we denote e| the 
dual edge oriented from the right face of Cy to the left face. The following notions 
are natural if we think of a dual 1-form on M as a 1-form on M* . A dual 1-form 
T is closed if for ever interior vertex i G V 

E^(4) = o- 

3 

It is exact if there exists Z : F —)■ M such that 

dZ{el) := Zyk - Zja = r(e;) 

where {ijk} denotes the left face of Cy and {jil} denotes the right face. 

We distinguish dual 1-forms from primal 1-forms for the following reasons. Firstly, 
the closedness conditions are different. The closedness conditions are imposed on 
faces for primal 1-forms while they are imposed at vertices for dual 1-forms. Sec¬ 
ondly, a discrete notion of the Hodge star operator is needed to identify 1-forms 
with dual 1-forms, although it is not explicitly used in this paper. In Discrete Ex¬ 
terior Calculus m one often uses the Hodge star operator, which maps a primal 
1-form w to a dual 1-form *a; via 

:= (cot -f cot 4)a;((e* )*) = -(cot -f cot I3'y)uj{ey) V{zj} S E 

where € (—tt, tt) denotes the angle Zjki of the triangle {ijk} with respect to 
some discrete metric, i.e. an assignment of edge lengths. Given a dual 1-form r 
and a primal 1-form df, we will occasionally write 

t(4) = kijdficy) 

for some k : Eint —t R. Here we think of it as r = fc' * df for some k' : E^t —t R- 

3. MOBIUS INVARIANCE 

In this section we prove that the class of isothermic triangulated surfaces is 
invariant under Mobius transformations. 

Given a triangulated surface / : F —>■ R^ and a Mobius transformation tr : 
R^U{oo}—^-R^Uloo}, we define cr o /: F —>■ as the triangulated surface with 
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vertices {a o f)^ := a o f^. We consider only the Mobius transformations that do 
not map any vertex to infinity. 

Taking a to be minus the inversion in the unit sphere, we obtain a triangulated 
surface 



Later we will identify with imaginary quaternions, which explains the notation 
f~^. We are going to show that / is isothermic if and only if f~^ is isothermic. 
We first rewrite the equations from Definition o 

Lemma 3.1. Given a non-degenerate realization f : V ^ of a triangulated 
surface, a -valued dual 1-form t : —>■ satisfies 

^T(e*) = 0 

3 

dficij) X r(e|) = 0 V{ij'} G Ei^t, 

Y^{df{eij),T{e*)) =0 Vz G 

3 

if and only if there exists k : Eint ^ such that 

kijdf{e^) = T(e*) V{zj} € Eint, 

^ ^ kjj df ( Cjj ) = 0 Vi € ^int: 

3 

3 

Proof. Suppose k : Ei^t ^ satisfies for every interior vertex i 

Yi hjdfie^j) = 0 . 

3 

Then, we have the identity 

^(d/(e^),M/(e^)) = - 2(/.- - /-/*)) = 

j 3 3 

Using this it is easy to verify all our claims. □ 

With the above lemma, we can show / is isothermic if and only if f~^ is isother¬ 
mic. 

Lemma 3.2. Suppose a non-degenerate realization f : V ^ of a triangulated 
surface is isothermic with a non-trivial dual 1-form r satisfying Definition \l . 1\ We 
write 

'^( 4 ) = kijdficy) 

for some k : Ei^t K. Then, the triangulated surface f~^ : U —>■ zs isothermic 

with corresponding dual 1-form 

G(e*) := d/”^(eg). 
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Proof. We check that f satisfies the equations in Definition o by applying the 
previous lemma. Firstly for every interior vertex i 

3 3 

- i/*i') + i/*i'E%(/* - fo) 

3 3 

= 0 . 

Secondly, for every interior vertex i 

E %i/*n/.f d/Ei' - i/r'n = E%(i/d^ - i^i') = 

3 3 

Hence, f~^ is isothermic with 1-form f satisfying Definition HU □ 

Proof of Theorem \l It follows from the previous lemma and the fact that Mobius 
transformations are generated by inversions and Euclidean transformations. □ 


Remark 3.3. The above calculation can be simplified if written in terms of quater¬ 
nions. Identifying the Euclidean 3-space with the space of purely imaginary quater¬ 
nions we obtain 


df H^ij) = fi ^df{eij)f^ 1 = /j. ^df{eij)f^ ^ 

and 

^( 4 ) = 

These two formulas are similar to the smooth case [33] . 

Lemma [3.I| provides another characterization of isothermic triangulated surfaces. 
We consider the light cone 

L ■.= {x €^3 \ x\ -\- x\ -\- x\ + x\ — xl = 0}. 

Corollary 3.4. Suppose f : V ^ is a non-degenerate realization of a tri¬ 
angulated surface and k : Ei^t —> K is a funetion. Then f is isothermic with 
corresponding dual 1-form t defined hy 

— kijdfi^Cij^ ^ Eijii 

if and only if for every interior vertex i 

(5) E kijdf{eij) =0 Vi e V^t 

3 


where / : H —>■ L C is the lift of f to defined by 


A function k : Emt —> K satisfying Equation ® is called a self-stress of /. 

It is known that the Mobius geometry of K^U {oo} is a subgeometry of the 
projective geometry of KP'^. Mobius transformations of M^U{oo} are represented as 
projective transformations of KP^ preserving the quadric defined by the light cone 
L. If two non-degenerate realizations are related by a projective transformation, 
then the spaces of self-stresses of the two realizations are isomorphic |53|. Hence, 


we obtain another proof of Theorem 1.3 
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4. Infinitesimal conformal deformations 

We consider infinitesimal conformal deformations for a given closed triangulated 
surface in space. We show that a surface is isothermic if and only if it is a singular 
point in the space of all surfaces conformally equivalent to the original one. 

4.1. Conformal equivalence of triangle meshes. We recall that a discrete met¬ 
ric of a triangulated surface is a function £ : E ^ K_|_ satisfying the triangle in¬ 
equality on every face. A non-degenerate realization f : V ^ induces a discrete 
metric £ : E ^ K+ via 


k'=\h-fi\ V{7j} e A. 

Definition 4.1 ([l^). Two discrete metrics £,£ : E - 
M are conformally equivalent if there exists u : V 

{y} 


1+ on a triangulated surface 
K such that for every edge 


£ij — e 


'£■■ 


Two non-degenerate realizations /, / : 14 —>■ are conformally equivalent if their 

induced discrete metrics are conformally equivalent. 

It leads naturally to an infinitesimal version of conformal deformations. 

Definition 4.2. An infinitesimal deformation of a non-degenerate triangulated 
surface f : V —>■ is a. map f : V —i' E^. It is conformal if there exists u : 14 —>• M 
such that the change of the induced discrete metric i : E ^ E satisfies for every 
edge {if} 

U, Ui „ 


£■■ — 




In particular, / is an infinitesimal isometric deformation if u = 0. 

The conformal equivalence class of a triangulated surface in Euclidean space is 
Mobius invariant [5]. It can be distinguished via logarithmic length cross ratios. 

Definition 4.3. Given a discrete metric £ : E ^ K+ on a triangulated surface, its 
logarithmic length cross ratio loglcr : —)■ is defined by 

loglcr(£)g ;= log^fe - log4i + log4 - log4 

where {ijk} is the left face of eij and {jil} is the right face. 

Theorem 4.4 ([5]). Two discrete metrics £ and £ on a triangulated surface are 
conformally equivalent if and only if 

loglcr(£) = loglcr(£). 

Corollary 4.5 ([5]). The dimension of the space of the conformal equivalence 
classes of a triangulated surface is \E\ — |I4|. 

4.2. Infinitesimal deformations. In this section, we consider closed triangulated 
surfaces. Suppose £ : E ^ K+ is a discrete metric on a closed triangulated surface. 
We consider an infinitesimal change of the discrete metric £ and write it as £ = a£ 
for some infinitesimal scaling a : E —>■ E. Then the change of logarithmic length 
cross ratio on edge {ij} is given by 

(loglcr(£))'- = ajk - cfki + cfu - crij =: L{a)y. 


The image of the linear map L : —>■ is the tangent space of the space of 

conformal equivalence classes (which is the same space at all discrete metrics). 
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Lemma 4.6. Given a closed triangulated surface. The operator L is skew adjoint 
with respect to the standard product {,) on KI'®! given by {a,b) := 
for any a,b G KI-®! . 

Proof. Let 5*^ : E ^Rhe the function defined by {S^^)ij = 1 on edge {ij} and zero 
on other edges. Then for any b G we have 

L*{b),, = {5\L*{b)) = {L{5^),b) = -bjk + bk^ - bu + by = -L{b),,. 

Thus we have L* = —L. □ 


The above lemma implies that we have an orthogonal decomposition 
— Ker(L) 0 Im(L*) = Ker(L) 0 Im(L). 

Lemma 4.7. Given a closed triangulated surface. We have the following. 

Ker(L) = {a : i? — >■ K | 3m G s.t. aij = Ui + Uj V{*j} G E} 

Im(L) = {a : i? —>■ K| ^ =0 \li gV} 

3 

Proof. It is obvious that 

{a : E ^ K|3m G s.t. Oij = Ui + Uj y{ij} G E} C Ker(L). 

Assume a G Ker(L). For each face we define 

_0 ajk 

(.Dj Ui := - 

Suppose Aiij is the neighboring triangle sharing the edge {ij} with A^k. Because 
of L{a)ij = 0 we have 

Uij 0 a^i Ujfc 0 an aij 

^ ^ = u,. 

Since the link of each vertex is a disk (although we only need the vertex link to be 
a fan), Equation ([^ in fact defines a function u : E —)■ M such that for any edge 
{*)■} 

dij — ■ 

Hence 

Ker(L) = {a : E —>• K | 3m G s.t. = Ui + Uj y{ij} G E}. 

On the other hand, it is obvious that 

Im(L) C {a : E —> K I Oy =0 Wi G V}. 


Since 

rank(L) = \E\ — dimKer(L) = \E\ — |E| 
the two vector spaces are indeed the same. 


□ 


Recall that conformal equivalence classes of a triangular mesh are parametrized 
by logarithmic length cross ratios. By the inverse function theorem the result below 
implies that by deforming a non-isothermic surface in space we can reach all nearby 
conformal equivalence classes. It is precisely in the case of an isothermic surface that 
the hypothesis of the inverse function theorem fails to be satisfied. Thus the space 
of all non-isothermic non-degenerate realizations in a fixed conformal equivalence 
class is a smooth manifold. 
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Theorem 4.8. Suppose f : V ^ is a non-degenerate realization of a closed 
triangulated surface. Then f is isothermic if and only if there exists a non-trivial 
element a € Im(L) such that 

(a,L(cr)) = 0 

for all infinitesimal scalings a : E ^ M. coming from infinitesimal extrinsic defor¬ 
mations in Euclidean space, i.e. for which there exists f '■ V —> : E ^ ISS' 

such that df — adf + df x W. 

Proof. Suppose / is isothermic with r satisfying Definition |l.l[ Let / : T —>■ be 
an arbitrary infinitesimal deformation and we write df = adf + df x W. Since r is 
closed, i.e. = 0 Vi S 1/ we have 

0 = - /i) = ^ydfidj) + df{ey) x Wy). 

iGV j {ij}eE {ij} 

From df{eij) x T(e*) = 0 we obtain 

0= XI +d/(ei,) X ^ (r(e*),fi/(eg))cr,j-. 

{ii}eE {ij}eE 

Using 

{T{e*),df{eij)) = (r(e*),d/(e^-*)) 

we see that (r, df) : E —>■ M. is well defined. Since we know that for every interior 
vertex i 

^(d/(e*j),r(e*)) = 0 
j 

we thus have {T,df) G Im(L). Hence there exists an non-trivial element a G Im(L) 
such that for every edge {ij} 

L{a)ij = -{T{e*j),df{ey)). 

Because / is arbitrary we conclude that 

0 = ((t, df),a) = (-L(a), a) = (a, L{a)) 

for all infinitesimal scaling cr : U —>■ K coming from infinitesimal extrinsic deforma¬ 
tions. 

On the other hand, suppose there exists a non-trivial a G Im(L) such that 

(a,L(cr)) = 0 

for all infinitesimal scaling a G coming from infinitesimal extrinsic deforma¬ 
tions. We define a dual 1-form r : E*^^ -G via 

df{eij) X T(e*) = 0, 

(d/(e,j),r(e*)) = -L{a)y 

for every edge {ij}. Since {df,T) G Im(L), we have 

X(d/(e^),T(e*))=0 WiGV. 

3 

In addition, for any infinitesimal deformation / : U —)■ we write df = adf-\-dfxW 

for some a \ E and W : E ^ R^. We obtain 

= (^’d.(a)) = 0. 

iev j ij ij 

Since / is arbitrary we conclude that r is closed, i.e. 

^r(e*) = 0 VfGU. 
j 
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Hence, / is isothermic with dual 1-form r. □ 

Proof of Theorem \1.4\ Consider the composition of maps 

{infinitesimal deformations in (infinitesimal scalings} (change of Icrs}. 

The space of infinitesimal conformal deformations is exactly Ker(L o cr). Moreover, 
we know 

dim(Ker(L o cr)) = 3|H| — rank(L o a) > 3|H| — (|i?| — |y|) = \V\ — 6g + 6. 

Finally we conclude: The inequality is strict L o cr is not surjective / 

is isothermic. □ 

Since the conformal equivalence classes of a triangle mesh are parametrized by 
length cross ratios, we can rephrase the previous theorems as follows. 

Corollary 4.9. Given a closed triangulated surface, isothermic realizations are 
precisely the points in the space of all non-degenerate realizations where the map 
that takes a non-degenerate realization to the conformal equivalence class of its 
induced metric fails to be a submersion. 

It is interesting to see how combinatorics affect geometry. It is known that the 
number of vertices of a closed genus-g triangulated surface satisfies the Heawood 
bound [18] 

|y| > 7 + Vl + 48g ^ 

This condition is known to be sufficient for the existence of a genus-g triangulated 
surfaces with \V\ vertices except for g = 2. Comparing the Heawood bound with 
the inequality in Theorem 1 1. 4| we obtain more examples of isothermic surfaces. 

Corollary 4.10. Every non-degenerate realization of a closed triangulated surface 
with \V\ < 6g -b 4 is isothermic. 

Proof. The space of infinitesimal conformal deformations contains all deformations 
that come from infinitesimal Mobius transformations. Therefore this space has 
dimension at least 10 and hence a surface must be isothermic if 10 > \V\ —6g-\-6. □ 

Some of these surfaces with small number of vertices can be realized in Euclidean 
space without self-intersection. For example, there are embedded surfaces with 
g = 2 and \V\ = 10 as shown in [^ . 

5. Preserving intersection angles 

Given a triangulated surface in Euclidean space, every triangle determines a 
circumscribed circle and every two triangles sharing an edge determines a circum¬ 
scribed sphere if the vertices are not con-circular. Two circumscribed circles are 
called neighboring if their corresponding triangles share an edge. We will call two 
circumscribed spheres neighboring if they have a common vertex. 

Intersection angles of circles and spheres are Mobius invariant. The intersection 
angles of neighboring circumcircles of a triangulated surface were used to define 
discrete Willmore functional |7]. 

Proof of Theorem \1.5\ Suppose we have an infinitesimal deformation / that pre¬ 
serves the angles between circumcircles and circumspheres but is not induced from 
Mobius transformations. Then it cannot be that / also preserves the length cross 
ratios (because it is not hard to see that in this case / is an infinitesimal Mobius 
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transformation). We write df = adf + df x W for some ct : —>■ K and W : E 

Then the change of logarithmic length cross ratios is L{a) where 

^jk ^ki T ^il ^ ddijii. 

(See Figure]^) By our assumptions L{a) does not vanish identically. 

We define a dual 1-form 

( ")■" ( )" M/(e,)p- 

Then we have 


'r(e*) X df{ey) = 0 y{ij} e E 

=0 Vi G V,nt- 

3 3 

In order to show that / is isothermic, we need to verify the closedness of r, i.e. for 
every interior vertex i 

51^(4) = 0- 

3 

We identify Euclidean space with the space ImH of imaginary quaternions 
(SectionjlT]). Pick any vertex vq and denote its neighboring vertices hy vi,V 2 , ■ ■ ■ ,Vn- 
Then we take an inversion in the unit sphere centered at /g ;= f{vo) and denote 
the images of the neighboring vertices by fi. We have the following relations: 

4 -f0 = if 3 - fo)-\ 

/, + ! - /, = -(/, - /o)-'((/, + l - /o) - ifj - foWj + l - /o)-' 

= -if3 - /o)-'(/l + l - /,){/, + ! - /o)-'. 


We define the infinitesimal scaling 


•- 


14+1 fj 


l/i+i-Z.I 

By taking the logarithmic derivative of the following equation 


14+1 ~41 ^ l/j+i ~/ill/j-i ~/ol 

l4 “ 4-il iZi+i “ /oll/j - /t-il 

we obtain for j = 1,..., n 


^3,3 + i- - ^3-^,3 - i^i^))0j 

where a^j = \ fj - fi\'/\fj - ft\- On the other hand, the vertices /i, 4, ■■■,/«, 4 
form a closed polygon in We define 

4j+i •= 14+1 “ 4l> 

4+1 - 4 


T, 


+1+1 


Since the polygon is closed, we have 


14+1-41 


^ — 54 4m+i4j+i- 
1=1 


The fact that the deformation / preserves the intersection angles of neighboring 
circles and neighboring spheres implies that the angles between the neighboring 
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segments and osculating planes of the closed polygon remain constant. Thus there 
exists a constant vector c G such that 

= ~ fj) 

= -^L(a)o,(/,-/o)-' 

=E-(4)- 

i=i 

To show that the converse is true one only has to reverse the previous argument. □ 

6. Example: Isothermic quadrilateral surfaces 

We show that isothermic quadrilateral surfaces as defined by Bobenko and Pinkall 
[2] are isothermic under our definition (after an arbitrary subdivision into trian¬ 
gles). Isothermic quadrilateral surfaces are analogous to conformal curvature line 
parametrizations of smooth isothermic surfaces. They can be treated using the the¬ 
ory of integrable systems. New isothermic surfaces can be obtained from a given 
isothermic surface via the Christoffel duality and Darboux transformations |20j . 
Special discrete surfaces related to isothermic quadrilateral meshes were studied in 

mm- 

Questions about infinitesimal rigidity of quadrilateral meshes have been consid¬ 
ered by [351 m] • 

We first review some results on isothermic quadrilateral surfaces from |S]. Then 
we construct an infinitesimal isometric deformation for every isothermic quadrilat¬ 
eral surface and show that the change of mean curvature around each vertex is zero. 
In this way we obtain isothermic triangulated surfaces from the earlier notion of 
isothermic quadrilateral surfaces. 

6.1. Review. 

Definition 6.1 (|2]). A discrete isothermic net is a map F : —>■ for which 

all elementary quadrilaterals have factorized real cross-ratios in the form 

Pn 

where am G R does not depend on n and /Sn G K not depend on m. 


Fm,n Fm-\-l^n 

Figure 1. An elementary quadrilateral 

Theorem 6.2 (0)- Let E : Z^ —^ fee a discrete isothermic net. Then the 

discrete net E* : Z^ —>■ defined (up to translation) by the equations 

^ Fm+l,n Fm^n 

-^m+l,n -^771,71 11 ^ ^ 112 ? 

||-^m+l,n 7n,n\\ 

P* _ F* _ R ^m,n+l ~ Pm,n 

m,n+l m,ra— Pniip '\\2 

\ |-^ m,n+l ^ 7n,n \ \ 
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is isothermic. F* is called the Christoffel dual of F. 

We need a formula for the diagonals of its Christoffel dual (Corollary 4.33 in [3). 


Lemma 6.3. Given a discrete isothermic net F, the diagonals of any elementary 
quadrilateral of its Christoffel dual are given by 




Pn)- 

M- 




F, 


m+l,n+l 


-Frr 


12 ’ 


Fm+l,n 

Ip _ p 112 ■ 

Dm+l.n 


6.2. Infinitesimal flexibility of isothermic quadrilateral surfaces. Given a 
discrete isothermic net we first arbitrarily introduce a diagonal for each quadrilat¬ 
eral in order to get a triangulation. Then, we define infinitesimal rotations on faces 
as follows. 

Rule: Suppose ABCD is an elementary quadrilateral of a discrete isothermic net 
F and the diagonal AC is inserted. Then we get two triangles ABC and 

ACD. We define infinitesimal rotations Zabc ■= B* and Zacd '■= D* where B* 
and D* are the corresponding vertices of the Christoffel dual F* : 


Theorem 6.4. Suppose we are given a discrete isothermic net F : Z^ — )■ and 
its Christoffel dual F* : 1? ^ MA . We assume that the faces of F have been 
subdivided into triangles in an arbitrary way. Then the infinitesimal rotations given 
by the above rule for each triangle define an infinitesimal isometric deformation 
of the triangulated surface. Moreover, the infinitesimal deformation preserves the 
integrated mean curvature and is not induced from Euclidean transformations. 


Proof. By Theorem |6.2| and Lemma |6.3[ the infinitesimal rotations of two adjacent 
triangles are compatible on the common edge. Therefore they define an infinitesimal 
isometric deformation. 

It remains to show that around an arbitrary vertex the change of the integrated 
mean curvature around vertices is zero. For every vertex there are 2^ = 16 ways of 
inserting diagonals on the four neighboring quadrilaterals. Taking the symmetry 
into account we can reduce them to 6 cases. We enumerate these 6 cases and 
calculate the change of mean curvature on each edge in Table It can be checked 
directly that in all cases the sum around the vertex is zero. □ 


Remark 6.5. Although the infinitesimal rotations on faces depend on the trian¬ 
gulation, the deformations of the edges already present in the quad mesh do not. 
For example, the change of the edge Fm+i,n — Fm,n is given by 




— (Frn-t-l.n F^n^n) ^ 


/T** _i_ P* 

■^m+l,n m^n 


Here we have used —|| {Fm+i,n — Fm,n)- Moreover, the quadrilaterals 

do not stay con-circular under the infinitesimal deformation. 


The infinitesimal isometric deformation defined above has an exact counter¬ 
part in the smooth theory. Given a simply connected isothermic surface / and 
its Christoffel dual f *, there exists an infinitesimal isometric deformation / satisfy¬ 
ing df — dfx f* (Section]^. It preserves the mean curvature but does not preserve 
curvature lines. If in addition the curvature lines were preserved, the shape operator 
would remain unchanged and the deformation would be trivial, i.e. an infinitesimal 
Euclidean transformation. 
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Table 1. The six types of triangulations around a vertex 
and the corresponding change of mean curvature on the edges. 


7. Example: Homogeneous discrete cylinders 

In this section we show that every homogeneous triangulation of a circular cylin¬ 
der in is isothermic. Here “homogeneous” means that there is a subgroup of 
Euclidean transformations that acts transitively at vertices and respects the com¬ 
binatorics. Note that in general none of the edges of such an isothermic discrete 
cylinder is aligned with the curvature line directions of the underlying smooth cylin¬ 
der. 

We consider the group G of all Euclidean motions that fix the z-axis. Every 
element g S G is of the form that acts on a point p as 

/ cos 9 sin0 0 \ f ^ \ 

9ip) = I — sin6l cos 9 0 ) P + ( 0 j 

Vo 0 1/ \hj 

where d, h € M. 

We pick two elements gi ,g 2 of G in general position and consider the group 
generated by gi,g 2 - For a generic choice of gi, (72 the group H is isomorphic to 1?. 
An element (s, t) G 1? corresponds to the element gfg^ G H. 
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We also consider 1? as the vertex set of a triangulated surface with faces of the 
form {(s, t), (s + 1, t), (s, t + 1)} or {(s + 1, t), (s + 1, t + 1), (s, t + 1)} (Figure]^. 



Figure 2. A triangulated surface with vertex set 1? 

We now define a map / : —>■ by picking r > 0 and setting 

/(s,t) = glgl{r, 0 , 0 ). 

For suitable 51,52 G G this map / will be a non-degenerate realization. Figure]^ 
shows a piece of such a discrete surface. 



Figure 3. A strip of an isothermic triangulated cylinder 


We now prove that realizations / : —>■ constructed as above are isother¬ 

mic. We will do this by showing that they admit a non-trivial infinitesimal iso¬ 
metric deformation preserving the integrated mean curvature. Note that up to 
symmetry there are only three types of edges, represented by {/(0,0),/(I, 0)}, 
{/(1,0),/(0,1)} and {/(0,1),/(0,0)}. We denote their lengths by 


4 (r, Oi, hi, 62,112) 


£b{r,6i,hi,92,h2) 


4(t, 01 , hi, 62 , h 2 ). 


The integrated mean curvature is the same at all vertices. We denote it by 


H{r, 01, hi, 02, h 2 ). 

Now the derivative of the map g := {£a,ib, 4, H) : K® —>■ has a non-trivial kernel 
at every point (r, 0i, hi, 02, ^ 2 ) G Moreover, it is easy to see that any non-zero 
element 

(r, 01, hi, 02, h 2 ) G kerd^ 

corresponds to an inhnitesimal deformation of / which is not induced from Eu¬ 
clidean transformations. This infinitesimal deformation preserves all the edge lengths 
and the integrated mean curvature around vertices. Therefore the triangulated 
cylinder / is isothermic. 
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8. Example: Planar triangular meshes 

In this section we show that certain planar triangular meshes are isothermic. 
For a simply connected surface, we know that a realization is isothermic if and 
only if there exists a non-trivial infinitesimal isometric deformation preserving the 
integrated mean curvature (Corollary |1.2[ ). We will see that every such deformation 
of a planar triangular mesh is given by a discrete harmonic function in the sense of 
the cotangent Laplacian [30l |37j . 

Theorem 8.1. Suppose f : V ^ C is a non-degenerate realization of a 
triangulated surface with Euler Characteristic x o-fid |I4| boundary vertices. Then 
f is isothermic i/ |Vb| — Sy > 0. 

Proof. Since each boundary component is a simple closed polygon, the number of 
boundary edges is \Eb\ = |I4|. The Euler characteristic is given by 

|E|-|ii;| + |F|=x. 

Since the surface is triangulated, we have 

3\F\ = 2\E\-\Eb\. 

Hence 


\E,nt\ - 3\V,nt\ = |14|-3x. 

By Lemma 3.1 a dual 1-form satisfying Definition o is equivalent to a function 
k : Eint —>■ such that for every interior vertex i 


(7) 

( 8 ) 


^ kijdficij) = 0 

3 


which is a system of linear equations. By simple counting and using the fact that 
/ is planar we obtain a lower bound for the dimension of the solution space 

dimjA: : Eint R satisfying 0(|8])} > \Eint\ - 2|Ei„t| - \Vint\ = \Vb\ - 3%. 

Hence / is isothermic if |t4| — 3y > 0. □ 


In particular the above theorem implies that every planar triangulated disk 
(y = 1) with more than 3 boundary vertices is isothermic. Since a disk is sim¬ 
ply connected, by Corollary |1.2| there exists a non-trivial infinitesimal isometric 
deformation preserving the integrated mean curvature. The following indicates 
how to obtain such infinitesimal deformations. 

Theorem 8.2. Let / : E —>■ R^ C R^ 6e a strongly non-degenerate realization 
of a triangulated surface with normal TV S and u : V —)■ R &e a function. 
Then the infinitesimal isometric deformation f := uN preserves the integrated 
mean curvature if and only if u is a discrete harmonic function in the sense of the 
cotangent Laplacian, i.e. for every interior vertex i 

E(cot/3| -h COt/3',)(Mj - Mi) = 0. 

3 

where G (—7r,7r) denotes the angle Zjki of the triangle {ijk} under the image of 
f (Figure\^. 
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Figure 4. Two neighboring and oriented triangles sharing edge {ij}. 


Proof. It is easy to see that for any function m : F —>■ K the infinitesimal deformation 
/ = uN preserves edge lengths. The infinitesimal rotation of a face {ijk} is given 
by 

T -f Uk^f 


^ijk — 


‘^^ij k 


where A^k '■= {{df{eij) x df{eik))/‘2,N) is the (signed) area of the triangle {ijk}. 
This follows from 

df(.^ij^ — — df X Zijk. 

To preserve the integrated mean curvature around vertices, the map Z : F ^ 
has to satisfy for every interior vertex i 


(9) 0— {df{eij), Zijk Zjii) — {df(ekj), Zijk). 

{ij}GE:i {ijk}GF:i 


Note that 


{df(,ekj).Z^,k) ={df{ekj). _ uMe,k) + n ^{e_k. ) + Ukdf{e,,)^ 

/ jj-/ 1 “ uf)df{eki) + {uk - Ui)df{eij) 

= {df[ejk), — -- -) 

= - cot Pij{Uj - uf) - cotP^iiuk - Ui). 

Thus ([^ is equivalent to saying that for every interior vertex i 


0 = ^ (cot + cot P\i)(Uj - Ui). 

{ij}CEA 


Hence the infinitesimal deformation uN preserves the integrated mean curvature if 
and only if u is a discrete harmonic function. □ 


In fact, every infinitesimal isometric deformation of a strongly non-degenerate 
planar triangular mesh is of the form u N for some function u modulo infinitesimal 
Euclidean motions. 

It can be checked that an infinitesimal normal deformation u N of a, planar mesh 
/ is an Euclidean motion if and only if u is a linear function, i.e. if there exists a 
vector a J- N and a constant c € M such that 

w = (a, /) + c. 


9. Example: Inscribed triangular meshes 


Since the notion of isothermic triangulated surfaces is Mobius invariant (Theorem 


1.3), our results for planar triangular meshes can be rephrased for triangular meshes 


inscribed in a sphere. Theorem |8.1| immediately implies the following. 
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Corollary 9.1. Suppose f \ V ^ is a non-degenerate realization of a triangu¬ 
lated surfaces with Euler Characteristic x l^fcl boundary vertices. Then f is 
isothermic i/ | 14 | — 3% > 0. 

In particular, it follows that every inscribed triangulated disk with more than 3 
boundary vertices is isothermic. This is analogous to the fact that disks immersed 
smoothly in a sphere are isothermic. 

On the other hand, to show that a triangulated surface is isothermic, we can 
look for a non-trivial infinitesimal isometric deformation preserving the integrated 
mean curvature. For an inscribed triangulated surface, we will see that it suffices 
to find an infinitesimal isometric deformation. 

The following lemma shows that for any triangular mesh inscribed in a sphere a 
dual 1-form t satisfying Equation ([^ and (|^ in Definition 1 1.1 1 will satisfy Equation 
Q automatically. 

Lemma 9.2. Given a non-degenerate realization f '■ V —>-8^ of a triangulated 
surface and a function k : Eint —t R, then for every interior vertex i 

^ kijdfieij) =0 kij\df{eij)\^ = 0 . 

j 3 

Proof. Since |/| = 1, we have 

= = 2(/„^%(/. - /,)) = o. 


Theorem 9.3. Suppose f : V ^ is a strongly non-degenerate realization of a 
triangulated surface M. Then every infinitesimal isometric deformation preserves 
the integrated mean curvature. 

Hence if f is infinitesimally flexible, then it is isothermic. If f is isothermic and 
M is simply connected, then f is infinitesimally flexible. 

Proof. Suppose an infinitesimal isometric deformation is given by a rotation vector 
field Z : F The compatibility condition implies that there exists k : Ei^t 

R such that on every interior edge {ij} 

{^ijk ^jil) ~ ^ijdf{eij') 

where {ijk} denotes the left face of cij and {jil} denotes the right face. The previous 
lemma yields for any vertex i € Vint, 

Hi = Y = 0 - 

3 

Hence the integrated mean curvature is preserved. □ 

Example 9.4. Jessen’s orthogonal icosahedron is obtained from a regular icosahe¬ 
dron by flipping 6 edges symmetrically without self intersection [211 [n]. Its vertices 
are exactly those of a regular icosahedron and hence lie on a sphere (Figure]^. It 
is known to be infinitesimally flexible and thus isothermic. 

Note that the property of being isothermic is Mobius invariant. 

Corollary 9.5. The infinitesimal rigidity of a non-degenerate simply connected 
triangulated surface inscribed in a sphere is Mobius invariant. 

We can regard a Mobius transformation of a triangulated surface inscribed in 
a sphere as being induced from a projective transformation of the ambient space. 
Then the above corollary is simply a special case of the projective invariance of 
infinitesimal rigidity [23] . 
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Figure 5. Jessen’s orthogonal icosahedron 
10. Discrete minimal surfaces 

In the smooth theory a minimal surface is the Christoffel dual of its GauB map. 
We give a definition of discrete minimal surfaces that is inspired by this fact. 

Definition 10.1. Given a non-degenerate realization / : F —>■ of a triangulated 

surface, a non-constant map /* : F —)■ is called a Christoffel dual of / if 

df{ey) X df*{e*) = 0 V{zj} G Ei^t, 

Y,{df{e^J),dr{el))=0 VzGF„t. 
j 

Since we know from the previous section that triangulated disks inscribed in a 
sphere with more than 3 boundary vertices are isothermic, the following definition 
is natural. 

Definition 10.2. Given a non-degenerate realization / : F —>■ of a triangulated 
surface, its Ghristoffel dual /* : F —)■ is called a discrete minimal surface with 
Gauss map /. 

Equivalently, a discrete minimal surface is a reciprocal-parallel mesh of an in¬ 
scribed triangulated surface. The combinatorics of a discrete minimal surface is 
that of the dual cell complex and each dual edge is parallel to the corresponding 
primal edge. Figure shows a discrete minimal surface together with its Gauss 
map. 

Starting with a planar triangulated disk / together with a discrete harmonic 
function u we obtain a dual 1-form r satisfying Definition 0 Using Lemma |3.2| 
we can apply a stereographic projection <f> to / in order to obtain a dual 1-form 
for $ o /. Then integrating the dual 1-form on the dual mesh yields a discrete 
minimal surface. 

Corollary 10.3. Every simply connected discrete minimal surface is given by a 
discrete harmonic function on a planar triangular mesh. 

As an example we take a triangulated square / in the (cc, yj-plane (Figure]^ left 
shows the stereographic projection of /). Notice that on a Delaunay triangulated 
disk D in the plane, a discrete harmonic function u is uniquely determined by its 
boundary values. The choice u\dD = xy leads to a discrete Enneper surface (Figure 
1^ right). As a second example, we consider a triangulated annulus centered at the 
origin with a cut along the positive x-axis. We solve for the discrete harmonic 
function u with boundary values either given by u\dM = log l^l or by u\dM = argz. 
We obtain a discrete helicoid and a discrete catenoid (Figure respectively. 
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Figure 6 . Right: A discrete Enneper surface corresponding to 
the discrete harmonic function u\dM = xy. The edges are parallel 
to the primal edges of the Gaufi image shown on the left. 



Figure 7. A discrete catenoid and a discrete helicoid. 


Theorem 110.31 indicates that the Weierstrass data needed to construct a discrete 
minimal surface consist of a discrete harmonic function on a planar triangular mesh. 
Its smooth counterpart is the Weierstrass representation for minimal surfaces, which 
states that each minimal surface locally is given by a pair of holomorphic functions. 
The choice of such a pair of holomorphic functions is equivalent to prescribing 
the Gauss map and the Hopf differential of a minimal surface. A similar formula 
for discrete minimal surfaces is elaborated in |28j . where a notion of holomorphic 
quadratic differentials is introduced to link discrete harmonic functions and discrete 
minimal surfaces by considering infinitesimal conformal deformations. Furthermore, 
it turns out that our definition of discrete minimal surfaces bridges the gap between 
earlier notions of discrete minimal surfaces m- 


11. Smooth analogues 


The main goal of this section is to prove the smooth analogue of Theorem 9.3 


Theorem 11.1. For every infinitesimal isometric deformation of an immersion 
f ■. M ^ , the mean curvature is preserved. 


Beyond this we also use the opportunity to review some known results on smooth 
isothermic surfaces that directly correspond to our discrete results. We rely on 
the treatment of smooth isothermic surfaces by means of quaternionic analysis as 
developed in [251 HU 132] ■ 
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Definition 11.2. Given two immersions / and / of a surface M in = Im H, / is a 
spin transformation of / if there exists a quaternion-valued function A : M —>■ H\{0} 
such that 

df = Xdf A. 

In this case, the normals N, N of /, / are related by N = X~^NX. 

Theorem 11.3. Given an immersion / : M —)■ = ImH and a function X : 
M —> 1H\{0}. Then Xdf X is a closed 1-form if and only if there exists p : M —)■ M 
such that 

df AdX = -pX\df\^. 

Suppose / : M —> is a smoothly immersed surface and X, JX G TpM form 
an orthonormal basis in principal directions at p G M. The corresponding principal 
curvatures ki , K 2 G R is given by 

dN{X,) = Kdf{Xi). 

The mean curvature Hp at p then satisfies 

df{X) X dN{JX) - dfiJX) X dN{X) = {m + K2)N = 2HN. 

Writing the above formula into quaternions yields 

df AdN = 2HN\dff. 

If / is a spin transformation of /, then the change of the metric and the mean 
curvature can be expressed as follows. 

Theorem 11.4. If f is a spin transformation of f given by X then the followings 
hold: 

(1) / and f are conformally equivalent since jd/p = |A|^|d/p. 

(2) We have d(XdfX) = 0 and hence 3p : M —)• M such that df AdX = —pA|d/P. 

(3) H\df\^ = {H + p)\X\^\df\^. 

Now we use the fact that every infinitesimal conformal deformation can be ex¬ 
pressed as an infinitesimal spin transformation. 

Theorem 11.5. Suppose / : M —)■ = ImH is a simply connected smooth 

surface and A : M —)■ H zs o function. Then there exists an infinitesimal conformal 
deformation / : M —>■ = Im H given by 

df = 2 Im(d/ A) 

if and only if there exists p : M —)■ M such that 

-df AdX = p\df\\ 

In particular 

(|d/n- =4Re(A)|d/|2, 

{H\dfn =ip + 2ReiX)H)\df\^. 

Definition 11.6. A smooth surface / : M —>■ = ImH is called isothermic if 

locally there exists a non-trivial ImH-valued closed 1-form t such that 

df At = 0. 

Theorem 11.7. A smooth surface / : M —)■ = ImH is isothermic if and only 
if locally there exists a non-trivial infinitesimal isometric deformation such that the 
mean curvature is unchanged. 
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Proof. Let / be isothermic and t be a 1-form satisfying Definition 11.6 The closed¬ 
ness of T implies that on any simply connected open set U C M there exists 
X : U ^ ImH such that 

T = dX. 


Since df A dX = 0, Theorem 11.5 implies that the 1-form lm{df X) is closed and 


hence there exists an infinitesimal deformation f : U ^ ImH satisfying 

df — 2Im(d/A). 


Because A is purely imaginary and df A dX = 0, from Theorem 


11.5 


we have 


i\dfn = o, 

{H\df\y=0 


and hence / is an infinitesimal isometric deformation preserving the mean curvature. 
The converse is proved similarly. □ 


The smooth analog of Theorem |1.3| is a classical result of isothermic surfaces 

Ha- 

Theorem 11.8. The class of isothermic surfaces is Mobius invariant. 

We need the global existence of r in order to relate it to the space of immersions. 

Definition 11.9. A smooth surface / : M —>■ = ImH is called strongly isother¬ 

mic if there exists a non-trivial Im H-valued closed 1-form t on M such that 

df At = 0. 

In addition, if t is exact and t = df* for some /* : M —>■ = ImH, then f* is 

called a Christoffel dual of /. 

The smooth analog of Corollary |4.9| in Section is known: 

Theorem 11.10 l|l()jl. Strongly isothermic immersions of a closed surface are 
the points in the space of immersions where the map from the space of immersions 
to Teichmuller space (which assigns to each immersion the conformal class of its 
induced metric) fails to be a submersion. 

Furthermore, Theorem |1.5| is analogous to the infinitesimal version of the follow¬ 
ing: 

Theorem 11.11 m)- Suppose f,f : M ^ S^ are two conformal immersions 
that share the same Hopf differential but do not differ by a Mobius transformation. 
Then f and f are isothermic surfaces. 

Finally, we establish a smooth counterpart of Theorem |9.3| in Sectionj^ We first 
need a lemma. 

Lemma 11.12. Let M be a surface with a Riemannian metric and / : M —>■ 
be an isometric immersion. Let X = g dfiY) -\- hN be an W-valued function on 
M where g,df{Y) and hN are its scalar, tangential and normal components. Then 
we can express df A dX in terms of standard operators from the vector calculus on 
M as 

—df AdX= [ — curlF -f df{Jgradg — AT -|- gradh) — ((divT) -|- 2hH)N] |d/p. 
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Proof. In the following we assume that X G TpM is an unit tangent vector. We 
first consider the scalar component. 

-df A dg[X, JX) = -df{X)dg{.JX) + df[JX)dg{X) 

= Ndf{dg{X)X + dg{JX)JX) 

= Ndf{gvd,dg). 

Then we consider the normal component. 

-df A d{hN){X, JX) = {i-df A dh)N - hdf A dN){X, JX) 

= Ndf {grad h)N - h{df {X)dN {J X) - df{JX)dN{X)) 

= df {grad h) — 2hHN. 

Finally we look at the tangential component. Notice that for an immersed surface 
in Euclidean space the induced Levi-Civita connection is given as follows: for any 
tangent vector field Y and tangent vector Z, 

dfiVzY) = d{df{Y)){Z) - {d{df{Y)){Z),N)N 
= d{df{Y)){Z) + {df{Y),df{AZ))N 
= d{df{Y))iZ) + {Y,AZ)N 

where A is the shape operator of the immersion /. We recall the definition of curl 
and divergent operator of a tangent vector field Y: 

div(y) : = {X,XxY) + {X,XxY), 
curl(y) : = {JX,VxY) - {X,VjxY) 

= -{X,XxJY)-{JX,VjxJY) 

= -div(jr). 

Collecting the above information we now obtain 
-dfAd{df{Y)){X, JX) 

= - df{X){df{XjxY) - {Y,AJX)N) +df{JX){df{XxY) - {Y,AX)N) 

= {X,VjxY) - {JX,VxY) - {JX,VjxY)N 

+ {-X, VxY)N - {AY, JX)df{JX) - {AY, X)df{X) 

= — curlE — (divy)A^ — df{AY). 


□ 


Proof of Theorem 11.1 Suppose f : M ^ is an immersion and A : M —>■ Im(IHI) 
induces an infinitesimal isometric deformation of /. Then writing A = df{Y) + hN, 
we have 

(iJ|d/|2)- = -dfAdX = - curlr + d/(-r +gradh) - {{divY + 2h))N. 
Comparing the imaginary part yields 

Y = grad h 


and thus 

{H\df\‘^)' = — curl(F) = — curl(gradh) = 0. 
Since (|d/P)' = 0, we have H = 0. 


□ 
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